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Abstract 

We consider the damped-driven KdV equation 

u - .... + n_ - ,uu. = .l.eS\ jud.^j,d.^,, 

where < < 1 and the random process r] is smooth in x and white 
in t. For any periodic function u(x) let / = (/i, I2, . . . ) be the vector, 
formed by the KdV integrals of motion, calculated for the potential 
u{x). We prove that if n(t, x) is a solution of the equation above, then 
for < t < z^-^ and ^ the vector I{t) = {Ii{u{t, ■)), h{u(t, ■)), ...) 
satisfies the (Whitham) averaged equation. 
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Introduction 

It is well known since the pioneer works of Novikov and Lax that the KdV 
equation 

U + Uxxx-QuUx^O, (0.1) 

defines an integrable infinite dimensional Hamiltonian system in a space 
of 27r-periodic Sobolev functions of order p > with zero meanvalue. It 
means that KdV has infinitely many integrals of motion Ii,l2, . . . , which are 
non-negative analytic functions on H^, and for any non-negative sequence 

1 — [Ii, I2, . . .) the set Tj = {u : Ij{u) = Ij Vj} is an analytic torus in 
HP of dimension |J(/)| < 00, where J is the set J = {j : Ij > 0}. Each 
torus Tj carries an analytic cyclic coordinate (p — {(pj,j G J{I)}, and in the 
coordinates (/, cp) the KdV-dynamics takes the integrable form 

7 = 0, ip = W{I) . (0.2) 

The frequency vector W analytically depends on /. See [MT76, KP03] and 
Section 2 below. 

Importance of these remarkable features of KdV is jeopardised by the fact 
that KdV arises in physics only as an approximation for 'real' equations, and 
it is still unclear up to what extend the intcgrability property persists in the 
'real' equations, or how it can be used to study them. 

The persistence problem turned out to be difficult, and the progress in 
its study is slow. In particular, it was established that small Hamiltonian 
perturbations of KdV do not destroy majority of time-quasiperiodic solu- 
tions, corresponding to (0.2) with |J(/)| < 00 (see [KukOO, KP03]), but 
it is unknown how these perturbations affect the almost-periodic solutions 
(I J(/)| = 00), and whether solutions of the perturbed equations are stable 
in the sense of Nekhoroshev. 

Probably it is even more important to understand the behaviour of solu- 
tions for KdV, perturbed by non- Hamiltonian terms (e.g., to understand how 
small dissipation affects the equation) . The first step here should be to study 
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how a t/-perturbation affects the dynamics (0.2) on time-intervals of order 
z/~^. For perturbations of finite-dimensional integrable systems this ques- 
tion is addressed by the classical averaging theory, originated by Laplace 
and Lagrange. During more than 200 years of its history this theory was 
much developed, and good understanding of the involved phenomena was 
achieves, e.g. see in [AKN89] . In particular, it is known that for a perturbed 
finite-dimensional integrable system 



where / e R", (/? e T", on time-intervals of order v ^ the action I{t) may be 
well approximated by solutions of the averaged equation 



provided that the initial data (/(O), </'(0)) are typical. This assertion is known 
as the averaging principle. 

The behaviour of solutions of infinite- dimensional systems (0.3) on time- 
intervals of order > is poorly understood. Still applied mathematicians 
believe that the averaging principle holds, and use (0.4) to study solutions of 
(0.3) with n — oo. In particular, if (0.3) is a perturbed KdV equation, writ- 
ten in the variables {I,^p), then (0.4) is often called the Whitham equation 
(corresponding to the perturbed KdV). The approximation for I{t) in (0.3) 
with < t < by /(t), satisfying (0.4), is called the Whitham averaging 
principle since in [Whi74] the averaging is systematically used in similar sit- 
uations. In so far the Whitham averaging for the perturbed KdV equation 
under periodic boundary conditions was not rigorously justified. Instead 
mathematicians, working in this field, either postulate the averaging prin- 
ciple and study the averaged equations (e.g., see [FFM80] and [DN89]), or 
postulate that the solution regularly - in certain sense - depends on the small 
parameter and show that this assumption implies the Whitham principle, see 



The main goal of this paper is to justify the Whitham averaging for 
randomly perturbed equations. 

Let us start with random perturbations of the integrable system (0.2) 
with I G M", (/9 G T", where n < oo. Introducing the fast time t — i/t 
write the perturbed system as the Ito equation 



I^uf{I,cp) 0^W{I) + iyg{I,cp), 



(0.3) 




(0.4) 



[KriSS]. 



dl ^ Fdr + adl3r, 

dip = {u-^W{I) +G)dT + g d(3r ■ 



(0.5) 
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Here F, G, a and g depend on (J, (/?), [3^- is a vector- valued Brownian motion 
and <T, g are matrices. It was claimed in [Kha68] ^ and proved in [FW03] 
that (under certain assumptions, where the main one is non-degeneracy of 
the diffusion a and of the frequency-map W) when z/ — > 0, the solution /(r) 
converges in distribution to a solution of the averaged equation 

dl ^ {F){I)dr+ {a){I)d/3, (0.6) 

where (F) is defined as in (0.4) and the matrix {(t){I) is a symmetric square 
root of the matrix J^„ aa^ d(fi. 

Now let us consider a randomly perturbed ('damped-driven') KdV equa- 
tion 

u - vu^x + u^xx - Quu^ = v^vit, ■ (0-7) 

As before, x E and J udx = J r]dx = 0. The force 77 is a Gaussian random 
field, white in time t: 

^ = ^ X] ^s(^s{t)es{x) , 

where Zq = Z \ {0}, (3s{t) are standard independent Wiener processes, and 
{csiS e Zo} is the usual trigonometric basis 

es{x) = { . ' (0.8) 

I smsx, s < . 

Concerning the real constants hg we assume that 

hs<Cm\s\~'^ ym,s (0.9) 

with some constants Cm (so i]{t,x) is smooth in x), and 

bs^Q Vs. (0.10) 

The factor y/u in front of the force 77 (i, x) is natural since under this scaling 
solutions of (0.7) remains of order 1 as t — * 00 and u ^ 0. Eq. (0.7) defines 
a Markov process in the function space H^. Due to (0.10) it has a unique 



""^The main theorem of [Kha68] deals with the situation when the unperturbed system 
is a stochastic equation with a non-degenerate diffusion for (p, but in its last section it is 
claimed that the ideas of the proof also apply to (0.5). 
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stationary measure. Let u'^it, x), t > 0, he a corresponding stationary in 
time solution for (0.7); or let u'^ be a solution, satisfying 

u^iO, x) =uo{x), (0.11) 

where Uo{x) is a non-random smooth function. In Section 1 we prove that all 
moments of all Sobolev norms ||w''(t, ■)\\m are bounded uniformly in i/ > 
and t > 0. Let us write u'^{r) as {I^{t), (p^{t)). These processes satisfy the 
infinite-dimensional equation (0.5), so by the just mentioned estimates the 
processes {/'^(•), < < 1} form a tight family, and along suitable sequences 
Uj ^ we have a weak convergence in distribution 

r^(-)-/0(-), (0.12) 

where, according to the type of the solutions u'^{t), the limiting process I^{t) 
is either stationary in r, or satisfies /°(0) = I{uq{-)). 

The main results of this work are the following two theorems, proved in 
Section 6: 

Theorem A. The limiting process /°(t) satisfies the Whitham equation 
(0.6), corresponding to the perturbed KdV equation (0.7). It is non-degene- 
rate in the sense that for any r > and each k > Iwe have P{I^{t) = 0} = 0. 

Theorem B. If the processes u'^{t) are stationary in r, then for any r > 
the law of the pair (/'^^ (r), ip'^^ (r)) converges to the product measure x dip, 
where g° is the law of /°(0) and dcp is the Haar measure on T°°. 

The proof is based on the scheme, suggested by Khasminskii in [Kha68] , 
see also [FW98] and [Vcr91]. It uses the estimates from Section 1 and more 
sophisticated estimates, obtained in Sections 4 and 5. Namely, we use cru- 
cially Lemma 4.3 (Section 4) and Lemma 5.2 (Section 5). In the former 
coupling arguments are evoked to prove that for any k probability of the 
event {Ik{t) < ^} goes to zero with S, uniformly in u and t. This is im- 
portant since (0.5) is an equation for I in the octant {J \Ij > Vj} which 
degenerates at the boundary {/ \Ij = for some j}. In the latter we examine 
the random process W"^{t) = !¥"*(/'' (r)), where W"^ is the vector, formed 
by the first m components of the frequency vector W. Exploiting Krylov's 
results from [KrySO] we estimate the density against the Lebesgue measure of 
the law of the averaged vector 14^™" (/'^(r)) dr, s ~ 1. We use this esti- 

mate to show that with probability close to one the components of the vector 
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W"^{r) are non-commensurable, so the fast motion (d/dr)(p'^ — i'~^W"^{t) 
is ergodic on the torus C T°°, for any m. This is a crucial step of the 
proof of Theorem A. Our proof of Lemma 5.2 is 'hard' in the sense that it 
uses heavily the analyticity of the frequency map W{I). 

The arguments above are applied to the perturbed KdV equation, written 
in the Birkhoff normal form (eq. (2.1) in Section 2). They apply as well to 
perturbations of other Birkhoff-integrable equations if their solutions satisfy 
good apriori estimates uniformly in the small parameter, and the correspond- 
ing transformation to the Birkhoff coordinates is smooth and is polynomially 
bounded at infinity. In the KdV case which we consider, half of the required 
bounds on the transformation is estabhshed in the recent paper [Kor]. We 
are certain that the remaining half can be obtained similarly, but do not 
prove them in this work, sec Theorem 2.3 in Section 2. 

The Whitham equation (0.6), corresponding to the perturbed KdV (0.7), 
is a complicated infinite-dimensional stochastic differential equation. Theo- 
rem A implies that for any smooth initial data 7(0) it has a weak solution, 
but we do not know if this solution is unique. We point out that, firstly, if 
(0.6) has a unique solution and the process 'u'^(r) satisfy (0.11), then the law 
of the limiting process /° is independent of the sequence {i^j}, and the con- 
vergence (0.12) holds for t/ — >• 0. Secondly, if (0.6) has a unique stationary 
measure, then a similar assertion holds for stationary solutions u'^{t). 

The inviscid limit. Let us consider the stationary solutions of eq. (0.7) 
in the original time t. The apriori estimates from Section 1 imply that this 
family is tight in C([0, T]; H^) for any T > and any p > 0. Therefore, along 
sequences — > 0, we have convergence in distribution 

(the limiting process u^{t) apriori depends of the sequence {I'j}). The argu- 
ments, applied in Section 10 of [KukOG] to the randomly perturbed Navier - 
Stokes equation (0.14) also apply to (0.7). They show that a.e. realisation of 
the limiting process u^{t, x) is a smooth solution of the KdV equation (0.1). 
In particular, the law ijP of the random variable m°(0, ■) G is an invariant 
measure for the dynamical system which KdV defines in H^. But KdV has 
infinitely many integrals of motion; so it has a lot of invariant measures. How 
to distinguish among them the measure /jPl Noting that u''{t)t=o — u'^{r)r=o, 
we apply Theorem B to get that the isomorphism u{-) i— > (7, </?) transforms 
to the measure g° x dip. In particular, if (0.6) has a unique stationary 
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measure, then the measure is uniquely defined, and the convergence (0.13) 
holds for z/ — > 0. 

This discussion shows that in difference with the deterministic situation, 
averaged randomly perturbed equations describe not only behaviour of solu- 
tions for a pre-limiting equation on time-intervals of order i/'^, but also its 
asymptotic in time properties. Indeed, under the double limit 'first i — > oo, 
next u 0', the distribution of any solution converges to a measure, simply 
expressed in terms of a stationary measure of the averaged equation. 

The Eulerian limit. The perturbed KdV equation (0.7) is a reasonable 
model for the randomly perturbed 2D NSE 

ii — vAu + {u • V)u + Vp = \/vri{t, x), x E T^, 

f f (0-14) 

divu — 0, udx= r]dx = 0, 

obtained by replacing in (0.14) the 2D Euler equation (0.14),^=o (which is a 
Hamiltonian PDE with infinitely many integrals of motion) by KdV. Under 
restrictions on the random force f]{t, x), similar to those imposed on the 
force in (0.7), eq. (0.14) (interpreted as a Markov process in the space of 
divergence-free vector fields u{x)), has a unique stationary measure, see in 
[Kuk06]. Let {u'^{t), V^if)) be the corresponding stationary solution. Then, 
along sequences i/j ^ 0, the convergence in distribution holds 

(A), /(•)), (0.15) 

where the limiting process (m°, p") is stationary in time, is sufficiently smooth 
in t and x, and a.e. its reahsation satisfies the free Euler equation (0.14),^=o- 
Accordingly, the law /i° of vP{Q) is an invariant measure for the dynamical 
system, which the Euler equation defines in the space of divergence-free vec- 
tor fields. To study the measure (in fact, the set of measures since it is 
possible that now the limit depends on the sequence {I'j}), is an important 
problem in (mathematical) 2D turbulence. The problem, addressed in this 
work, may be considered as its model. 

Agreements. Analyticity of maps Bi — > B2 between Banach spaces Bi and 
B2, which are the real parts of complex spaces BI and Sg, is understood 
in the sense of Prechet. All analytic maps which we consider possess the 
following additional property: for any R a map analytically extends to a 
complex > 0)-neighbourhood of the ball {|m|bi < R} in B^. When two 
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random variables are equal almost sure, we usually drop the specification 
"a.s.". 

Notations, xa stands for the indicator function of a set A (equal 1 in A 
and equal outside A) . By x{t) we denote various functions of t such that 
x(t) — > when i — > oo, and by Xoo{t) denote functions x{t) such that 
>c{t) = o{t~^) for each N. We write x(t) = x(i; R) to indicate that x(t) 
depends on a parameter R. For a measurable set Q C we denote by \Q\ 
its Lebesgue measure. 



1 The equation and its solutions 

We denote by H the Hilbert space 

H^{ue L2{S^) : J udx^O} 

with the scalar product {u,v) = ^ J^^ u{x)v{x) dx. Then {es,s G Zq} (see 
(0.8)) is its Hilbert basis. We set H"^ to be the m-th Sobolev space, formed by 
functions with zero mean- value, and given the norm = (f^; i^)^^^- 
We write the KdV equation as 

u + V{u)^0, V{u) ^ u^^^ - 6uux , (1.1) 

and re- write eq. (0.7) as 

u - vu^x + V{u) = ^/iyr]{t, x) . (1.2) 

It is well known that a dissipativc nonlinear equation in one space- 
dimension with a white in time r.h.s. has a unique strong solution if the 
equation's solutions satisfy sufficiently strong a-priori estimates. In Appendix 
we show that any smooth solution of (0.7) with a deterministic initial data 

u{0)^uo, (1.3) 
where Uq G H"^,m > 1, satisfies the following estimates: 

Ee-ll«(*)llo < max(Ee"ll"(°)ll'>, 2e2'^^o), (1.4) 

E\\u{t)\\l<ma^{4E\\um\l,C'J, (1.5) 

E\\u{t)\t < C{\\uo\\mk,B^+^,m,k) . (1.6) 
Here i > 0, A; G N and a < (2 maxfe^)-!. 
Accordingly, we have the following result: 
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Theorem 1.1. For any deterministic Uq G H'^, m > 1, the problem (0.7); 
(1.3) has a unique solution u{t,x). It satisfies estimates (1.4) - (1.6). 

Due to assumption (0.10), eq. (0.7) has a unique stationary measure Hi, 
and any solution converges to /Xj, in distribution. For the randomly forced 
2D NSE equation this result now is well known (e.g., see in [Kuk06]). The 
proofs for eq. (0.7) are simpler and wc do not discuss them. 

Let M°(t, x) be a solution of (0.7), (1.3) with uq = 0. Since V{u1{t)) Hu, 
then Theorem 1.1 and the Fatou lemma imply 

Theorem 1.2. The unique stationary measure iiy satisfies the estimates 
I e'^ll"llo/x^(rfM) < < oo V(T < (2max62)-i, 

JH 

Jh 



2 Preliminaries on the KdV equation 

In this section we discuss intcgrability of the KdV equation (1.1). 

For r > let us denote by h^' an abstract Hilbert space with the basis 
{fj,j = ±1, ±2, . . . } and the norm | • 1^, where 

\v\l-j:r''{v] + vl,) for v^J2-jfr 



We denote v 




, and identify a vector v = Yli'^jfj ^ with the 



sequence (vi, V2, . . . ). 

Theorem 2.1. (see [KP03]). There exist an analytic diffeomorphism ^ : 
H ^ hP and an analytic functional K on of the form 

K{J2 ^3fj) = K{h, I,,...), /, = i {v^ + vl,) , 

with the following properties: 

1) ^ defines, for any m an analytic diffeomorphism : if"* — > hJ^, 

2) d*(0) is the map H"^ 3 ^u^Ss ^ E \s\~^^^vjs G h"^; 
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3) a curve u{t) e (7^(0,7; if) is a solution of (1.1) if and only if v{t) — 
'^{u{t)) satisfies the equations 

Vj = -sign(i) /2, . . . ) , i e , (2.1) 

where Wi ^ for I ^ 1,2, ... . 

Corollary 2.2. If u{t) is a solution of (1.1) and ^{u) — v — J^'^sfs, then 
hit) ^l{vl + v^,){t)^ const yk^l,2,.... (2.2) 
li V & h^, then the vector / = (Ji, /2, . . . ) belongs to the space 
h^j^{I: |7|,. = 2^/+2'-|/,.|<oo}. 
In fact, I e h}_^_, where 

h'i+^{lehy. Ij>0 Vj}. 

Amplification. The function K in Theorem 2.1 is analytic in That is, 
it analytically extends to the vicinity on this set in the space h^. 

The quantities /i,/2, . . . are called the actions. Each vector Vj can be 
characterised by the action Ij and the angle 

ipj = arctan — - . 

We will write v — {I, (f), where (f — {(pi, (f2, ■ ■ ■)■ The vector (f — {(fi, (f2, . . . ) 
belongs to the torus T°°. Wc provide the latter with the Tikhonov topology, 
so it becomes a compact set. 

The functions u — >■ Vk{u), k e Zq, form a coordinate system on H. They 
are called the Birkhoff coordinates, and the system of equations (2.1) - the 
Birkhoff normal form for the KdV equation. The normal forms is a classical 
tool to study finite-dimensional Hamiltonian systems and their perturbations 
locally in the vicinity of an equilibrium (see [MS71], §30). For all important 
finite-dimensional systems the normal forms do not exist globally. In con- 
trast. Theorem 2.1 shows that the KdV equation is an infinite-dimensional 
Hamiltonian system which admits a normal form globally in the whole space 
H. To take all advantages of this normal form we will need some information 
about asymptotic properties of the transformation ^(m) when u ^ oo: 
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Theorem 2.3. For m = 0, 1, . . . there are polynomials Pm and Qm such that 

W^{u)U<P^{\\u\U), J = 0,1,2, 

and 

\\d^'^'\v)\\m<Qm{\v\m). J = 0, 1, 

for all u,v and all m > 0. Here for j > I |<i-'^'|m is the norm of the 
corresponding poly-linear map from H"^ to K^, and similar with ^d^'^~^^ra- 

Proof. The estimates for the norms and follows from 

Theorem 2.1 in [Kor].^ 

We do not prove here the estimate for d^'^{u) with j — 1,2. We are certain 
that modern spectral techniques (e.g., see [Kor, DM06]) allow to establish 
them, but wc think that this paper is not a proper place for a corresponding 
rather technical research. □ 

Remark. We do not use that the coordinate system v = (vi, V2, . . . ) is sym- 
plectic, but only that it puts the KdV equation to the form (2.1). Therefore 
we may replace v by another smooth coordinate system v' = (v^, V2, . . . ) 
such that Ij — Ij for all j and cpj = ipj + $j(/i, /2, ■ • ■ )■ Non-symplectic 
coordinate systems are easier to construct, and it is possible that a proof of 
Theorem 2.3 simplifies if we replace there f by a suitable system v'. □ 

For a function / on a Hilbert space H we write / e LipLock(-f^) if 

\f{u^)-fiu2)\<PiR)\\u,-U2\\ if \M, \\U2\\<R, (2.3) 

where P is a continuous function (depending on /) . Clearly the set of func- 
tions LipLock(-f^) is an algebra. Due to the Cauchy inequality any analytic 
function on H belongs to LipLock(-f^) (see Agreements). In particular, 

W^/ e LipLock(/i^) for ZeN, r>0. (2.4) 

3 Equation (0.7) in the BirkhofF coordinates 

For A; = 1, 2, . . . we denote 

^Notc that the quantity, denoted there || J||p_i/25 equals l^lp+i up to a constant factor, 
and Q2p satisfies the estimates Q2p < -Rip(||u||p+i) and ||w||p+i < R2p{Q2p), where Rip 
and R2p are some polynomials. 
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where ^(w) — v — (vi,V2, ...). Let u{t) = n^{t) be a solution of (0.7), 
which cither is a stationary solution, or satisfies (1.3) with a //-independent 
non-random uq. Applying Ito's formula to the map we get: 

Let us denote 

d*fe(«)(^6,e,d/3^) =5fc(«)ci/3=^5it,■(«)(i/3^ 5^,- e VA;,j. 

j 

Then the diffusion term in (3.1) may be written as y/u Bk{u) df3. 

Since = | is an integral of motion (see (2.2)), then application of 
Ito's formula to the functional | jv^p = and eq. (3.1) results in 

dh = y({d^!k{u)u:cx,^k) + ^ (X]&Jc?^^fc(w)[ej,ej], Vfc) 

1 X ' ^3.2) 

+ 2 &|M^fc(w)eirjc?t + x/^(Efc(M) d(3, Vfc) 

(here and below (■, ■) indicates the scalar product in M^). Note that in differ- 
ence with (3.1), eq. (3.2) 'depends only on the slow time' in the sense that 
all terms in its r.h.s. have a factor u or s/v. 

Let us consider the infinite-dimensional Ito process with components 
(3.2), k>\. The corresponding diffusion is y^crd/?, where cr = [akj{u),k G 
N,j e Zo) and 

CTkj = {Bkj{u),Vk) = bj{d'^k{u)ej,'^k{u)). 
Consider the diffusion matrix a, 

a{u) = a{u)a\u) , au^k.^ = '^kijcrk2j ■ (3-3) 

Lemma 3.1. For any u & H the sums in (3.3) converge. The matrix a is 
symmetric and defines a bounded linear operator in P . If — for some 
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^ G P, then ^fe 7^ only if = 0, where v = In particular, if 

Vfc 7^ V/c, then Ker a = {0}. Moreover, if \vj\ > S for 1 < j < m, then for 
any ^ e x {0} C M°° we have 

{a{u)^,Oi.-W'ml>C\^\l, (3.4) 
where C depends on 5,m, \v\i and the sequence {bj}. 

Proof. Using (0.9) and Theorem 1.1 we get that \crkj\ < C\j\^^rik, where 
7] & P. Therefore a defines a bounded hnear operator H ^ and ci* defines 
a bounded operator ^ H. So a = aa* is a bounded operator in and its 
matrix is well defined. Let us take any vector ^. Then (a^,^);2 = {a^i,a^i) 
where 

i^'Oj = bj{d^kiu)ej, ^rk)ik = bj{ej, d<i/{uy{®^k^k)). (3.5) 

k 

Hence, ^ eKera if and only if rf^(M)*(©CA;Vfc) = 0. Since d'i/{u) is an iso- 
morphism, then in this case ^^^k = for each k, and the assertion follows. 

To prove (3.4) we abbreviate ®^feVfe = ^„ and denote d^{u)*^v — rj. Then 
(7*(ii)^ —dia,g{bj}ri (see (3.5)). Due to the first assertion of Theorem 2.3, 

ll^ll?<C^l(l^|l)le.|?<C^l(|^|l)|e|?,C^2(|^|o)m^ 

So 

00 

k=N+l 

for any A^. Since {d^{u)*)~^ — (o?^(w)~^)*, then the second assertion of the 
theorem implies that 

00 

= ll^llo ^ ^o(klo) \U > Co{\v\o) lulls'. 

k=i 

Choosing iV= [(2CiC2C-ir2m3)'/']+lwegetthat Zk=iVl > |<^o(^'ICI/V 
Accordingly, 

k=l 

where C" depends on the sequence {bj} and N. □ 
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Wc see that the infinite-dimensional Ito process (3.2)fcgN, defined for I e 
has non- degenerate diffusion outside the boundary dh^_^ — U '■ — ^ 
for some j > 0}. 

By applying Ito's formula to the k-th angle (pk = arctan ( ) {k > 1) 
and using (2.1) we obtain 

d(Pk = Wk{I) + iy\vk\'^{d'^k{u)u^x: ) 

oo 

j=i jeZo 
+V^\^^k\-\Bkiu),vi)dp, 

where = ^ ^ . Denote for brevity the drift and diffusion coefficients 

in the above equation by Wk{I)+iyGk{v) and y/i' gl{v) respectively. Denoting 
similarly the drift coefficients in (3.2) by i'Fk{v) we rewrite the equation for 
the pair {Ik, (pk) {k > 1) as 

dlk{t) = uFk{v)dt + y/v(Tk{v) dPt , 

dcpkit) = [Wk{I) + iyGk{v)]dt + V^Qkiv) dPf 



Introducing the fast time 

T — vt 

we rewrite the system (3.6) as 



dlkir) = Fk{v)dT + ak{v) d^r , 

rl 1 (3-7) 

dMr)^ [-Wk{I) + Gk(v)\dr + gk(v)d(3r. 

Here P = {Pj,j £ ^o), where /3j(r) are new standard independent Wiener 
processes. 

In the lemma below Pk and PkN are some polynomials. 

Lemma 3.2. For A; e N, j e Zq we have: 

i) the function Fk is analytic in each space , r > 2 (so Fk E LipLock(^'^) ), 
and has a polynomial growth as \v\k — > oo; 
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ii) the function (7kj{v) is analytic in h"^ , r > 0, and for any N > 1 satisfies 

\'Jkjiv)\<J-''PkN{\v\r)yveh'; 

ii) for any r > 2,6 > and N > 1 the functions Gk{v)x{i^>s} o-nd 
gkj{v)x{ik>s} fl^e bounded, respectively, by 6'^Pk{\v\r) and 6^^j^'^ PkN{\v\r)- 

Proof. The assertions concerning the functions Fk and Gk follow from The- 
orem 2.3 since the set of analytical functions with polynomial growth at 
infinity is an algebra. To get the assertions about cxfc and Qk we also use 
(0.9). □ 



4 More estimates 

In this Section and in the following Sections 5-6 we consider solutions of 
equation (3.6), written in the form (3.7), which either are stationary in time, 
or satisfy the i/-independent initial condition (1.3), where for simplicity uq is 
smooth and non-random, 

m 

First we derive for these solutions additional estimates, uniform in u. 

Lemma 4.1. For any u > 0,T > and m,N eN the process I{t) satisfies 
the estimate 

E sup |/(r)|;^y. = E sup \v{t)\^^ <C{N,m,T). (4.1) 
o<T<r o<T<r 

Proof. For the sake of definiteness we consider a stationary solution v{t) = 
{v^(t)}. Cauchy problem (3.7), (1.3) can be considered in the same way. 
Applying Ito's formula to the expression k^I^ gives 

d{k-I^) = k-{{NI^-^Fk{v) + \n{N - l)/f-^ Y.^BkAv\^kf)dT+ 

+NI^-^ak{v)d(5r). 

Therefore, 

E sup k'^Ikir) < EA;"*/f (0) + 
o<T<r 
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r 

„ ^ CO 

+k"^E sup / (iV/f + -N{N - J2 <) 

0<T<T J ^ ^ ^ ' 
n J=l 



+ 



/t"E sup 

0<r<T 



r 

y iV/f-i(s)afe(^;)(i/33| < C(m,iV,T) . 



Doob's inequality, Lemma 3.2 and Theorem 1.2 have been used here. This 
relation yields the desired estimate. Indeed, by the Holder inequality we get 

AT 



sup |/(r)|^^)=2^E sup (y.\if'-^%{T)) 

{oo oo 

oo 

<CjvE sup f Vj'^('™+'^/f(r)') <Ci(m,7V,r). 



jV"--' 



o<T<r . , 

j=i 



□ 



In the further analysis we systematically use the fact that the functionals 
Fk{I,ip) depend weakly on the tails of vectors (p — {(pi, (p2, . . .). Now we 
state the corresponding auxiliary results. 

Let / e LipLock(^"0 ^ud v e /i", n > ni. Denoting by 11^, M > 1, the 
projection 

\j\<M 

we have \v — IlMv\ni < M~^'^~'^'^^\u\n- Accordingly, 

\f{v) - f{UM{v))\ < P{\v\n)M'^--^^l (4.2) 

Similar inequalities hold for functions on k], and (2.4) with r — implies 
that 

\Wk{I) - W,{UmI)\ < Pfe(|/U)M-". (4.3) 

The torus acts on the space 11^//^° by linear transformations ^Om^ 
9m e T^, where sends a point vm = {Im,^m) to {Im,^m + Om)- 
Similar, the torus T°° acts on hP by hnear transformations $e : {!■, ^) ^ 
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{I,(p + 9). The transformation ^0 continuously depends on ^ e T°°, in the 
strong operator topology. 

For a function / G LipLock(^"0 ^ define the average of / in 

the first N angles as the function 

(here id stands for the identity transformation in the space h*^ IIjv/i^), and 
define the average in all angles as 

{f){v)= [ f{^ev)de, 

where d6 is the Haar measure on T°°. The estimate (4.2) readily implies that 
|(/)iv(^) - (/)(^)| < P(it:)7V-("-"^) if \v\n <R. (4.4) 

Let V — (/, ip). Then {f)N is a function, independent of ipi, . . . , ipN, and 
(/) is independent of cp. I.e., (/) can be written as a function (/)(/). 

Lemma 4.2. Let f G LipLock(/i"^)- Then 

i) The functions {f)N{v) and {f){v) satisfy (2.3) with the same polyno- 
mial as f and take the same value at the origin. 

a) They are smooth (analytic) if f is. Moreover, if f is smooth, then 
(/)(/) is a smooth functions of the vector (/i, . . . , Im) for any M. If f{v) is 
analytic in the space h"'^, then (/)(/) is analytic in the space hj^. 

Proof, i) Is obvious. 

ii) The first assertion is obvious. To prove the last two consider the 

function g{ri, ra, . . . ) = (/)(vi, V2, . . . ), v^- = (^^^ ^ . Then g{r) = (/)(/), 

where = |rf for each I. The function g is smooth and even in each r^, 
j > 1. Any function of finitely many arguments with this property is known 
to be a smooth function of the squared arguments, so the second assertion 
holds. 

Now let f{v) be analytic. Denote by f}"^ the space of all sequences r = 
(ri, r2, . . . ) such that the corresponding vector v belongs to and provide 
it with the natural norm. If f{v) is analytic, then {f){v) also is analytic and 
g{r) extends analytically to an even function in a complex neighbourhood O 
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of i)^^ in C. This neighbourhood may be chosen to be invariant with 
respect to all involutions 

(ri, r2, . . . , r^, . . . ) (n, rs, . . . , -rj, ...), j = 1,2,.... 

The image Oj of O under the map 

(ri,r2,...) ^ (H'H'---) 
is a neighbourhood of h]^ in the complex space h^^ C. The function 

g{±^/2h,±V^2,...) =:9{Vl) 

is a well defined locally bounded function on Oj.''' For any its restriction 
to Of = Oi n nAr(/ij^ ®C) is a single- valued algebraic function on a domain 
in C^; so g{\/l) is analytic on Of for each A^. Hence, g{\/l) is analytic 
on Oi (see Lemma A.4 in [KP03]). Since g{Vl) = (/)(-?"), then the result 
follows. □ 

Let {I^{t), ^"{t)) be a solution of (3.7). In the lemma below we show that 
the processes /^(t), k > 1, do not asymptotically approach zero as — > 
(concerning the notation >i{6~^; M,T), used there, see Notations): 

Lemma 4.3. For any M eN and T > we have 

P{mi^4-(T) <S}< k{5-';M,T), (4.5) 

k<M 

uniformly in u > and < r < T. 

Here the difficulty is that the scalar process /^(t) = ||vfc(r)p satisfies 
equation (3.7), where the diffusion ak degenerates when vanishes. The 
equation for the vector-process Vfe(r) (see (7.1) below) has a non-degenerate 
diffusion, but its drift has a component of order u'^. To prove the lemma's 
assertion wc construct a new process Vfe(r) such that |vfc(r)| = |vfc(r)| and 
Vfc satisfies an Ito equation with a nondegenerate diffusion and coefficients, 
bounded uniformly in u. Then Ik = ||vA,.(r)p meets estimate (4.5) by a 
Krylov's theorem. The problem to perform this scheme is that the process 

^i.e., it is bounded uniformly on bounded subsets of O/. 
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Vk is constructed as a solution of an additional diffusion equation which is ill 
defined when Vk vanishes. We cannot show that the event 



{vi,{t) = for some < r < T} 

has zero probabihty and resolve this new difficulty by means of some addi- 
tional (rather involved) construction. 
For a complete proof see Section 7. 



5 Averaging along Kronecker flows. 

The flow 

where W e M°°, is called a Kronecker flow. In this section we study averages 
of functions f{v) — f{I, (f) along such flows. That is, we study the quantities 



1 

^ J f{I,ip + W"^t)dt, T>0. 



Lemma 5.1. Let f e LipLock(^"0> = i^^^) & h^-,n > ni > 0, and f is 
analytic in the space /i"^ . Then for each R' > 0, m & N and 5 > there is a 
Borel set n'^,{S) C {x G : |a;| < R'} such that < S, and for any 

^ n'^,{S), IVV^I < R' the estimate 



T 



i I f{I,^ + W^t)dt-{f)iv] 



< ^co(m, R', \v\n, f) + m-("-"i)p(|^;U) , 



holds uniformly in ip G T°°. Here P is the continuous function from (2.3) 
and is identified with the vector {W"^, 0, . . . ) e R°°. 

Proof. Let us first assume that f{v) = fijlmv) (i.e., v depends only on 
finitely-many variables). Then / = /(J™, 99'") is analytic in (p'^ and the 
radius of analyticity is independent of /, satisfying |7|;jO < R'. Now the 
estimate with P := is a classical result (e.g., sec in [MS71]). In general 
case we write / as / o 11^, + (/ — / o 11^) and use (4.4). □ 
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We will apply this lemma with W"^ = W"^{I), where I — I{t) is the I- 
component of a solution of (3.6). To do this we have to estimate probabilities 
of the events {W^™(/(r)) G 01^, {5)}. To state the corresponding result we 
introduce more notations. For any events Q and O we denote 

PQ(e») = p(Cgnc»), 

and 

Eq(/) = E((1-xq)/). 
Abusing language, we call Pq a probability. We fix any 

P>1, 

denote 

Br^{I: < it:}, 

and for > 1 consider the event 

= { sup \v''{t)\p>R], 

0<T<T 

where v'^{t) is a solution. Noting that < R' = R'{R,m) outside the 

event Qr, we denote 

fl{5) ^fl^,{5), R' ^R'{R), 0<5 <1. 

Finally, for M > m and < 7 < 1 we define 

= {/ e /ir , : min L < 7} . 

' ^ l<j<M ■' ^ 

Lemma 5.2. There exists M = M{R,m) > m such that 
j\n^{{W'^{I{s)) e m}W{s) e Q,}) ds < H{d-';R,m,^,T) , (5.1) 

uniformly in u > 0. ^ 

"^We recall that >c{t; R, m, 7, T) stands for a function of t which goes to zero when 
f — > 00, and depends on the parameters R, m, 7 and T. 
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Proof. Consider the function D{I) — det (^dWp/dlr ■ I < j,r < m). It is 
analytic in (see Amplification to Theorem 2.1), and D ^ since -D(O) = 
C"^, C ^ (see [KukOO], Lemma 3.3, and [KP03]). For a finite non- decreasing 
sequence of natural numbers a = {a^ < ■ ■ ■ < a^) we denote 



\a\ — a 



.N 



[a] = N 



and define the derivative d°'D{I)/dI°' in the natural way. 

Step 1: Study of the sets {/ e : \D{I)\ <e}, < £ < 1. 
By the analyticity any point /' e Br has a neighbourhood O C /i° such 
that 



97" 



> c V/ e O, 



where the sequence a — {a^ < • • • < a^) and c > depend only on the 
neighbourhood. Since Br is a compact subset of h% we can cover it by a finite 
system of neighbourhoods Oj, j — 1, . . . ,L, as above, where L — L{R, m). 
Then 



{I^Br: 
Let us denote 



«1, j = l,...,L} = 0. 



(5.2) 



M = max \aj\ , A?" = max [aj\ 



i<j<i 



and consider the sequence 

£ = £o < £i < • • • < Sat < 1 , £j 
where < £ < 1. Note that 



Sj e]l^ = £^2 ) for < J < iV . 



For m < [aj] we set 

21- = {/ e i3« : 



d 



d 



-D{I) <£„}. 



In particular, 21° . = 21° = {/ G iS^ : \D{I)\ < e] for each j. 
For < £ <^ 1 relation (5.2) implies that 

L 

21° = U ((21° \ 2ti^.) U (2li^ \ 2l^J U • • • U (2lL7l-^ \ 21L7I)) . 
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Step 2: An estimate for the integral Pq^{|D(/(s))| < e} ds. 
Due to the last displayed formula, the integral to be estimated is bounded 
by a finite sum of the terms 

'"^Pn«{/(«)e2t^,\2i;+^}(i«, r<K]. (5.3) 



/ 

Jo 



To estimate (5.3), we abbreviate ^ . . . -^D{I) — f{I). Then 



\ -{I^Bn: |/(/)| < and > • (5.4) 



Consider the Ito process z{t) = /(/(r)). We define the Markov moment 
r' = min{T > : |/(t)|/jP > i?^} A T, and re-define z{t) for r > r' as a 
continuous process, satisfying 

dz{T) = dl3l for T > r' . 

Since t' > T outside Qr, then outside VIr we have ;2(t) = /(/(r)) for < 
T <T. For 2;(r) we have 

dz{T)^c{T) dr + ^6j(r) d/?^ , 

where |c| < C{R,m), bj — 5j^i for r > r' and 6j = J2§i^'^kj for "T" < 
Denoting a = ^6^, we have a = '^{cr(J*)jk^jfVkf- So |a(T)| < C{R,m). 
Prom other hand, (3.4) in Lemma 3.1 implies that 

M 

|a(T)|>C(i?,m,7)^(V,/)2 if /(r) ^ . (5.5) 

i=i 

Applying Theorem 2.3.3 from [KrySO] to the process z{t), we get 

E / {X{\z{T)\<er}\(^{r)\dT <C{R,m,T)er. 
Jo 

By (5.4) and (5.5) the integrand is > £^_^iC(i?, m, 7) iil{t) e {%.\K^^)\Qi- 
Hence, 

rPnAHs) e {W,^\%;')\Q,}ds<ere;^,C{R,m,j,T) 
Jo 

= £(2-^)C(i?,m,7,r). 
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We have seen that 



/ PnMD{I{s))\<e}\{I{s)eQ,})ds<e^''''^C,{R,mn,T). (5.6) 

Step 3: Proof of (5.1). 

We have an inclusion of events 

{W"^{s) e m} \ IHs) e Q,} c [[{W"^{s) e m} \ iiHs) e Q,} 
U{D{I{s)) < e})) U {{\D{I{s)) < 8} \ {I{s) e Q,}) . 

Probabihty of the second event in the r.h.s. is ah^cady estimated. To estimate 
probabihty of the first event we apply the Krylov estimate to the process 
W"^{s). Re-defining it after the moment r' (see Step 2) and arguing as when 
deriving (5.6) we get that 



fT 



<\n{6)\'/"^C{R,m,j,e,T). 



Pnn {{W^'is) e m} \ iUis) e Q,} U {D{I{s)) < £})) ds 



Finally, choosing first e so small that the r.h.s. of (5.6) is < e and next 
choosing 5 so small that the r.h.s. of (5.7) is < s, we see that the l.h.s. of 
(5.1) is < 2e for any £ > 0, if 5 is sufficiently small. □ 

6 The limiting dynamics. 

Let us fix any T > 0, an integer p > 3 and abbreviate 

hP = h, hPj = hi, hPj^ = hi+ \I\hP = \I\ , \v\p = \v\ . 

Due to Lemma 4.1 and the equation, satisfied by I^ij), the laws C{I^{-)} 
form a tight family of Borel measures on the space C([0, T]; /i/+). Let us 
denote by Q° any its weak limiting point: 

Q'=\^u.C{F^{.)}. (6.1) 
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Our aim is to show that Q° is a solution to the martingale problem in the 
space hi with the drift operator {F){I) — (F2)(/), . . . ) and the 

covariance {A){I) = {{Aki){I)}, where 

j 

By Lemmas 3.2 and 4.2 the averages (Fj) and (^fc^) are analytic functions 
on hi- The covariance (A) is non-degenerate outside the boundary of the 
domain h^_^ in the following sense: let ^ e M-^ C M°° and / e |/| < i?. 
Then 

E > ^Kl/' if l^.l > 7 > for i < M, (6.2) 

where C > depends on M, R and 7. Indeed, the estimate follows from (3.4) 
with V = (/, ip) by averaging in ip. 

Our study of the limit Q° uses the scheme, suggested by R. Khasminskii 
in [Kha68] and is heavily based on the estimates for solutions v'^{t), obtained 
above. 

First we show that for any k the difference 

T 

4(t) - J {Fk){I{s))ds (6.3) 


is a martingale with respect to Q° and the natural filtration of cr-algebras. 
A crucial step of the proof is to establish that 



:= E max 

0<T<T 



{F,{F{s),cp''{s))-{Fk){F{s)))ds 







(6.4) 



as 1/ — > 0. Proof of (6.4) occupies most of this section. 
Let us fix an integer 

m > 1, 

denote the first m components of vectors I'^ and cp'^ by I"''^ and (p"''^, and 
rewrite the first 2m equations of the system (3.7) as follows 

dF'"" = F'^{r,ip'')dT + a'^(r,ip'')dPr, 

= (^-pF'"(r) + G""(r,(^'^)j(ir + 5f™(r,(^'^)(i/3^. 
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Here and afterwards we identify the vectors {I^, 1^,0,0, .. .) with I^''^, 
and the vectors {(p'(, . . . , (p'^^, 0, 0, ... ) with (f^''^. 

Denote = {Fk)m{I,<f )i={im,o),ip=o- By Lemma 3.2 there is a 

constant Ck{R) such that for any v = {I,(p), \v\ < R, we have 



\FkiI,^)-Fk{I"',if^)\<Ck{R)m-\ 

\{Fk)Un-{F'k){I)\<Ck{R)m-'. 
Define the event as in Section 5. Due to Lemma 4.1 



(6.6) 
(6.7) 



(here and in similar situations below the function x is //-independent). Since 
by Lemma 3.2 the function F^ has a polynomial growth in v, then this esti- 
mate implies that 

IE max [ Ffc(v^(s))(is-En„ max f Ffe(v^(s)) dsl < Xoo(-R). 

0<T<T Jq V V // ^0<T<T Jq V V // I V / 

The functions Fkil"'"" , ip"'"^) , (Ffc)^(/'^'™) and (Ffc)(/'^'™) satisfy similar rela- 
tions. So we have 



< Xoo(i?)+Ea« max 
o<r<r 



{F,{r{s), ^''{s))ds - Fk{r^^{s), ^''''^{s))}ds 



-FEq max 
o<T<r 



-FEn„ max 

" 0<r<T 



{F,(r'-(s),(^^'-(s)) - (F,)^(r'™(s))}(is 



mu{r^-{s))-{Fu){r{s))]ds 



< 



< x^{R)+Ck{R)m-^ 

-FEq^ max 



0<T<T 



T 

y" {F,(r'-(5),^'^'™(5)) - {Fk)^{r^"\s))}d, 
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The last inequality here follows from (6.6)-(6.7). It remains to estimate the 
quantity 



max 

0<T<T 



(1 - Xn^) ■ 



To do this we consider a partition of the interval [0,T] to subintervals of 
length uL, L > 1 by the points 

Tj = vto + vjL, < j < + 1 , 

where tk+i is the last point Tj in [0, T\. The constant L such that 



L>2, L<-v-^ 
- 2 



(6.8) 



and the (deterministic) initial point to G [0, Ij) will be chosen later. Note 
that 

-T < K -vL <T. 
2 - 



Denote 



J (^Fk{r'"'{s),ip'''"'{s))ds-{F,U{r^"'{s)))ds, 0<1<K. 
n 

Since outside the event Q,r we have 



F,(r'-(s), vp'^'-(s)) - (Ffe)^(r'-(s)) )ds 

for any r' < r" such that r" — r' < uL, then 



< iyLC{R) 



En„ max 

0<r<T 



y (F,(r'-(s),(^'^'-(s)) - {F,)m{r'"'{s)))ds 





(6.9) 
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To calculate the contribution from the integral over an ^-th subinterval, we 
pass there to the slow time t — y~^T. Now the system (6.5) reads as 



rf<^-'™(t) = (iy™(r) + z/G™(r, ^'')) dt + V^g'^ir, (p'')df3t . 

Denoting tj = Tj/u = to + jL we have: 



(6.10) 



+1^ 



<^ y{Ffe(r'-(x),(^^'"^(x))- 
ti 

-Fk [r^"\ti), ip^'^iti) + w"\r{ti)){x - tl)) ]dx 
{Ffe(r'-(iO, ¥^^'"^(^0 + w^{r{ti)){x - tl)) - {F,)^{r''-{ti))}dx 



{{F,)UF^"^{ti)) - {Fk)m{I'''"'{x))}dx 



To estimate the integrals T^^ — Tf we first optimise the choice of to- Defin- 
ing the event fl{5), the number M{R,m) and the set as in Section 5, we 
have 

Lemma 6.1. The non-random number to e [0, i^L) (depending on v and 6) 
can be chosen in such a way that 



1 ^ 

5Z P^^i < Xoo(i?) + x(7"'; R, m) + x(r^ 7, R, m) 

1=0 



(6.11) 



for all < S,^ < 1, where 

£i = ^rU {I{ri) e Q,} u {W^'iri) e n{6)} . 
Proof. Due to Lemmas 5.2 and 4.3, 

£ p(nRU {I{t) e Q^} u {W"\t) e n{6)}) dr 

< Xoo(-R) + x(7~^; R, m) + >c{8~^\ R, m, 7) . 
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Writing the l.h.s. as Yli=o^{^i) d'to, where £i is defined in terms of 
Ti — to + i^jL, and applying the meanvalue theorem we get the assertion. □ 

Applying the Doob inequahty and Lemmas 3.2, 4.1 to (3.6) we get that 

Pn«( sup \r{t)-r{ti)\>P{R)iyL + A) 
ti<t<ti+i 

< P( sup iy\ f a{y{s)d^^ > A^) < Cjv(i^L)^A-'^ , 
ti<t<ti+i Jti 

for all N and A. Choosing in this inequality A = (i/L)^/^, using (6.8) and 
denoting 

Qi^{ sup \r{t)-r{ti)\>P,{R){iyLy/'}, 

ti<t<ti+i 

where Pi is a suitable polynomial, we have 

Pa«(QO<>^oo((^^^)-';m). (6.12) 

Let us set 

Then (6.11) implies the estimate 
1 ^ 

j^J2p^i< >^oo{R) + x(7-^; R, m) + x(5-^; 7, R, m) 
1=0 

+ Koo{{i'L)~'^\m) =: K. 
Since Fk{I, (f) has a polynomial growth in I, then 

K I ^ 

^|(E-E^jT^|<P(i?)- J]P^,<«: (i = 1,2,3), (6.13) 
1=0 1=0 

where we denoted by k, another function of the same form as above. So it 
remains to estimate the expectations Ejr^T;^ and their sums in /. 
First we study increments of the process ip^^'^it). Let us denote 

^-,rn^t) - v^'^iM) - W^{F{ti)){t - ti) =: $r(i), ti<t< ti+1. 
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Then 

+ / g'^d(3,=: J1 + J2 + J3. 
Jti 

Outside the event Ti the term Ji estimates as follows 

|Ji| < P{R,m){vLf'^L. 
To estimate J2 and J3 we assume that 

P{R){vLfl' . (6.14) 

Then outside we have 

\m)\>\l Vie[i^,t,+i], k<m, 
so by Lemma 3.2 and (6.14) there we have 

I J2I < uLC{R)-f-^ < C\R){uLf/^. 
To bound J3 we introduce the stopping time 

t' = min{i > ti : min /^(i) < 7 or \F{t)\ > R} ^ U+i . 



Then 

rt'At 



pt'At 

{1 - XT,)\J3{t)\ < / r(s)d/3,|=: J^(i). 
We have i/E J^* Igf'^pcis < vL^~^C{R^m). So the Doob inequality implies 

P^<{ sup I J3I > (^^)'/'} < P{ sup 141 > (z/L)^/-^} 
ti<t<ti+i ti<t<ti+i_ 

< {vLf'^^-^C{R,m). 

We have seen that 

P^,{*r > ^'(^, m)i/^/3L^/3} < (i/L)V3^-iC(i?, m). (6.15) 
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Now we may estimate the terms T^. 

Terms T^. Since G LipLock(^), then by (6.15) 'probabihty' Pjf, that 
the integrand in T^^ is > C(-R, m)z/^/^L^/^ is bounded by (i/L)^/^7~^C(i?, m). 
Since outside J-'i the integrand is < C{R,m), then 



Terms T^. By Lemma 5.1, outside Ti 

Tf < v5-^C{R,m) + Lvm-^C{R). 

So 

^E^,Tf < {5L)~^C{R,m)+m-^C{R). 
I 

Terms T^. By Lemma 4.2, outside J^i we have Tf < P{R){vLf/^{vL). 

So 

J]E^jf <P(i?)(z.L)V3. 

Now (6.13) and the obtained estimates on the terms T| imply that 

^ E|77z| < K + i/^/^C(i?, m, L, 7) + ((5L)-iC(i?, m) + m-^C{R). 
I 

Using (6.9) we arrive at the final estimate: 
Sl'^ < x^{R)+C{R)m-^ + i/LC(i?) 

(6.16) 

+ (same terms as in the r.h.s. above) . 

It is easy to see that for any e > we can choose our parameters in the 
following order 

so that (6.8), (6.14) hold and the r.h.s. of (6.16) is < e. 
Thus, we have proved 

Proposition 6.2. The limit relation (6.4) holds true. 
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In the same way one can show that 



E max 

0<i<T 



{F,{r{s),^-'{s))-{F,){r{s))}ds 



asi/^0. (6.17) 



I From Proposition 6.2 taking into account the a priori estimates we 
finally derive 

Proposition 6.3. The process (6.3) is a square integrahle martingale with 
respect to the limit measure Q° and the natural filtration of a-algebras in 

c{[o,ooy,hi+). 

Proof. Let us consider the processes 

(r) = (r) - r{F,){n{s))ds, re[0,T], j = l,2,.... 
Jo 

Due to (3.7) and (6.4) we can write N^^ as 

Here M^^ — I^J — J Fk{I''\ (p'^^) is a martingale, and S^-' is a process such 
that 

E sup |S^-' (r)|^0 as i/j ^ . 

0<T<T 

This convergence implies that 

hm £(7V;^ (-)) = hm £(M,^^ (-)) (6.18) 

in the sense that if one limit exists, then another one exists as well and the 
two are equal. 

Due to (6.1) and the Skorokhod theorem, we can find random processes 
J^^ir) and J(t), 0<t<T, such that CJ^^i-) = Cr^{-), CJ{-) = Q°, and 

^ J in C([0,T],/ij) as i/^ ^ , (6.19) 

almost surely. By Lemma 4.1, 

P{ sup |r^(r)| >R}< CR-^ 

0<T<T 
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uniformly in Vj. Since (F^) G Lipj^ockC^/) by Lemmas 3.2 and 4.2, then (6.19) 
implies that the left limit in (6.18) exists and equals (6.3). By Lemmas 3.2 
and 4.1 the family of martingales M^^(r) is uniformly integrable. Since they 
converge in distribution to the process (6.3), then the latter is a martingale 
as well. □ 

Denote Zk{t) = Ik{t) — (Ffc)(/(s))(is. Using the same arguments as 

t 

above and (6.17) we can show that Zk{t)Zj{t) — J {Akj){I{s))ds is a Q°- 



martingale in C{[0,T);hi+). Combining the above statement we arrive at 
the following theorem, where T > and p > 3 are any fixed numbers. 

Theorem 6.4. Let the process ti^(i), < u < 1, be a solution of equation 
(1.2) which either is stationary in time, or satisfies the u -independent ini- 
tial condition (1.3), where Uq is non-random and smooth. Let ^(■u'^(r)) = 
v'^{t) = {I''' {t) , ip" {t)) . Then any limiting point of the family C{I^{-)} 
as V ^ Q is a measure in C{0,T; h^_^) which satisfies the estimates 

[ sup \I{r)\f^m Qo{dI{-)) <C{N,m,T) <oo VA^,meN 

J 0<T<T ^ 

and solves the martingale problem in C{0,T;h^) with the drift {F){I) and 
covariance {A) [I). 

Let o"°(/) be a symmetric square root of {A){I) so that {a^{I))a^{lY = 
{A){I). We recall that {A) (I) is a positive compact operator for each I e h. 

Corollary 6.5. Any limiting measure Q° as in Theorem 6.4 is the distribu- 
tion of a solution I{t) of the following stochastic differential equation 

dl = (F(/)) dr + a\l)dWr , (6.20) 

where Wt is a cylindrical Brownian motion on h'j. 

Proof Denote by f) the Hilbert space of sequences {xi,X2 ■ ■ ■ , } with the norm 

00 

= It is easy to check that f) is continuously embedded in 

/ij, thus all the coefficients (F(/)), a°(/) and {A)(I) are well-defined for any 

I el). 

By Theorem 6.4 and Lemma 4.1 the measure is concentrated on 
C{0,T;hf+*). Since this space is continuously embedded in C(0, T; ()), then 
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Q° is also concentrated on C(0, T;f}). Therefore, Q° is a solution of the 
above limit martingale problem in the Hilbert space i). It remains to use 
Theorem IV.3.5 in [Yor74] (also see [DZ92]). □ 



The limiting measure Q° and the process I{r) inherit the uniform in v 
estimates on the processes I'^{r), obtained in Sections 1-4. For example, 

P{4(r) <6} < >i{6-^] k) uniformly in r e [0, T], for any k>l. (6.21) 

In particular, /(r) G \ dh^j^ a.c., for any r > 0. 

Remark. Equation (6.20) is the Whitham equation for the damped-driven 
KdV equation (0.7). Our results show that it has a weak solution in the space 
for a given /(O) which is a deterministic vector in the space hf^ = nh^j_^_. 
In fact, the same arguments apply when /(O) is a random variable in hj^ 
such that E||7(0)||^ < oo, where N and p are large enough. □ 

Now we assume that u'^{t) is a stationary solution of (1.2). Then the 
limiting process /(r) as in (6.20) is stationary in r. We denote q'^ — £(7(0)) 
(this is a measure on the space hFj_^_). 

Theorem 6.6. Let a process u'^{t) be a stationary solution of equations (1.2). 



1) for any < r < T the law of (fi^ir) converges weakly as u ^ to the 
Haar measure dip on T°°. 

2) The law of the pair {I^{r), ip^{T)) converges, along a subsequence {vj}, 
corresponding to the measure Q^, to the product measure q'^ x d(f. 

3) For any m the measure (f'^ = C{I"^{0)) is absolutely continuous with 
respect to the Lebesgue measure on ]R!{! . 

More precisely, the second assertion of the theorem means the following. 

Due to (6.21) the limiting measure us supported by the Borel set (1 {v : 
Vj 7^ Vj}, which is measurably isomorphic to (/ij^ \ dh^_^_) x T°°. Under 
this isomorphism the limiting measure reeds as g° x dip. 

Proof. 1) Let us fix any m and take a bounded Lipschitz function /, defined 
on the torus c T°°. Then 



Then 




T 



1 



E / /(^^'-(5))d8, 



33 



where ip^'"^ satisfies (6.5). Arguing as when estimating the expectation in 
the Lh.s. of (6.9) in the proof of Theorem 6.4, we get that 

Jo 

Therefore E/((/7'^(t)) —>(/), and the first assertion of the theorem follows. 

2) Consider an arbitrary bounded Lipschitz test function of the form 
$(/,(^) = /(/™)y((^'"), m > 1. We have 

E$(r(r),9.^(r)) = ^^f ^ f{r^"'{t))g{ip-^'"^it))dt. 

Consider a uniform partition of the interval (0, i/~^T) into sufficiently long 
subintervals. As was shown in the proof of Theorem 6.4, with high probability 
on any subinterval of the partition the function /'^'™(t) does not deviate 
much from a random constant (see (6.12)), while the normalised integral of 
g{(p'^'"^{t)) approaches the integral of g against the Haar measure (see the 
proof of the first assertion). Therefore when i/ — > 0, the r.h.s. above can be 
written as 

(^E J f{F'-'{s))ds) j gi^^d^ + o{l) = 

T°° 

= j fin dq' j g{^^) d^ + o{l) . 

This completes the proof of 2). 

3) The vector I™-{t) satisfies the Ito equation, given by the first m 
components of (6.20). The corresponding diffusion is non-degenerate by 
(6.2). Therefore by the Krylov theorem (see [KrySO]) for any Borel set 
U C [5, S-^Y", 5 > 0, we have that 

g°"^(C/) = ^{F'^it) eU}< Q|C/|^/'". (6.22) 

Let us take any zero-set Z C M™ and write it as 

Z = Zi U • • • U U Z , where Zj C {Ij = 0} and Z C M^g . 

Then q^'^iZj) = for each j due to (6.21). Writing Z = U5>o^5, where 
Z5 = Zn [5,5-^]™, we use (6.22) to get that g°"^(Z) = Yirsiq^'^{Zs) = 0. So 
(f'^[Z) — and the theorem's proof is completed. □ 
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Remcirk. For any j > 1 the measure g° = £(/j(0)) satisfies an analogy 
of estimate (6.22) with m = 1. Therefore g° = fj{s)ds, s > 0, where the 
function fj is bounded on segments [5, 

7 Proof of Lemma 4.3 

Step 1: processes v^(t). 

For rii,ri2 G \ {0} we denote by U{rji,rj2) the element of 5*0(2) such 
that [/(7yi,r;2)^ = ^. Note that U{r]2,r]r) = t/(r/i, r^a)"' = [/(r^i, 7^2)*- 

In the fast time r equation (3.1) reads 

dVfc= (id*fe(ii)y(ii)+Afc(^;))(iT + ^Sfe,(^;)(i/3^ (7.1) 

j 

where we denoted 

Ak{v) = d^!k{u)u^:c + ^ X] h]d'^^!k{u)[ej, Cj] , Bkj{v) = d^!k{u)hjej . 

Let v{t) = {vji.(r), A; > 1} be a solution of the system (7.1)^^^. 
We introduce the functions 

Ak{vk:v) = C/(vfe, Vfe)Afe('y), Bkjiyk.v) = U{vk:Vk)Bkj{v) , 

smooth in (vfe, v^) from (M^ \ {0}) x (M^ \ {0}), and consider the additional 
stochastic equation for Vfe(T) e M^: 

dvk = Ak{vk, v)dT + Bkj{vk, v)d(5l . (7.2) 

j 

Its coefficients arc well defined for all non-zero and v^. 

If v{t) = {vfc(r), /c > 1} is as above, then eq. (7.2) with a prescribed 
initial data has a unique solution, defined while 

|vfe|, |vfe| > c, < C, 

where c, C are any fixed positive constants. This solution may be obtained as 
the last component of a solution {v, v^) of the coupled system (7.1)^^^! (7-2). 
This system has a unique solution since (7.1)^^^ is equivalent to (0.7) (so it 
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has a unique solution), while (7.2) is a Lipschitz equation on the domain, 

defined by the conditions above. 

For a 7 e (0, |) we introduce the stopping times , i > and r^", i > 1, 
where Tq = Q and for i > 1 

T-=inf{T>T+i : |vfc(T)| <7 0r \v{t)\^, > , 

T+ = inf{T > 7^- : |vfe(T)| > 2^ and \v{t)\^, < i-)} . 

Note that Tq < rf , t[' < < t^j^-^ if i > 0, and t"^ — >• oo as j — cxd. 

Next we construct a continuous process v^(t), t > 0. We set ^],{tq) = 
Vfe(ro"). For i = 1 we extend v^(t) to the segment Aj_i := [t^-i^t^] as a 
solution of equation (7.2), and on the segment Aj = Vi' I'^'t] define it as ^ 

v^(t) = C/(v,(7v-), Vfc(T-))vfc(T), for T e A, . (7.3) 

Lemma 7.1. -(/" |v^(Tj^i)| = |vfc(r^"^]^)| and vl satisfies (7.2) on Aj_i, t/ien 
|v^| = |vfc| everywhere on that segment. 

Proof. Application of Ito's formula to the expression II = \ |v^P on the 
segment Aj_i yields 

dil = (v^, ^(v^, v))dT + 5] (1 \B,,{^rl v) \'dT + (v^, S,,(v^, v))d(3\) . 

I 

Similarly, Ik — ||vjkp satisfies 

dh = {wk,Ak{v))dT + Y^ \Bki{v)\^dT + {wk,Bki{v))di3[) . 

I 

By construction, the drift and diffusion coefficients of these two equations 
satisfy the relations 



^If v/j(0) = 0, then Tq'" =0 and the formula (7.3) is not defined. But it happens 

with zero probability, and in this case we simply set = 0. 
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For the squared difference {1^ — we have 

d{h - ilf = (2(4 - p J^'\~f'\ ^k,A,{v))+ 



1 (Iv.l 










2 



X (7-4) 
"^{wk,Bki{v)f) dr + dMr, 
I 

where Air is a square integrable stochastic integral whose structure is of no 
interest. Denote JT(r) = (4 - /fe)^((T V t+) A r^i)). Since 



7-7 

, , , ^ ^k - J 



T — T"i 
-71 _ o -'fc 



I I I ~ 7 1 ' 



then it follows from (7.4) that EJT(r) < EJt(O) + C(7) /EJt(s)c^s. As 







J"'{r^i) = 0, then J^{t) = by the Gronwall lemma. That is, |v^| = |vfe| 
on Aj_i. □ 

Applying this lemma with i = 1 we see that (7.3) with i = 1 is well 
defined, and |v^| = |vfc| on Aq U Ai. Repeating the construction above for 
i — 2,3, . . . we get a continuous process v^(t), r > 0, satisfying (7.2) on the 
segments > 0, satisfying (7.3) on the segments Aj,i > 1, and such that 

|v^(r)|^|v,(r)|. 

Let us abbreviate f/j = L'"(vfc(rj~), Vfc(r~)). Then on the intervals Aj the 
process v^(t) satisfies the equation 

dvlir) = Ui(^^d^k{u)V{v) + Ak^dr + UiBkj{v)d/3i. 
Finally, using the notation 

i 

U,^ld^k{u)V{v)+Ak) r e U(Tr,T+), 



Ak{^rk,v,t) 



and 

{Bkj{vk,v), relJAi, 
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we represent v^(t) as the Ito process 

T r 

^lir)=Vk{0) + J A,{^lv,s)ds + J BkA^lv,s)d(3l. (7.5) 



|v^l Iv'l 

Letting formally — 1 for |vjt| = 0, we make the function = 1 along 
all trajectories. 

By the definition of and B^j and by Theorem 2.3 the following bounds 
hold true with a suitable integer K: 



\Ak\<C{\v\'l+l), reUA, 
|i,|<C^-niHf + 1), ^eU(7;-,r+) 

i 

\Bk\h^ < C(|^;|f + 1), re [0,cx)) 

(cf. Lemma 3.2). Let us fix any z/ > 0. The family of processes {v^(-) , < 
7 < 1/2} is tight in C(0,T;M^). This readily follows from (7.5), Lemma 4.1 
and the estimates above. 

Since Bkj{v) — /3jd^fe(w)ej, where ^ defines diffeomorphisms — > hP 
and ^ h^, then the diffusion '^Bkjdp^ in is non-degenerate and the 
corresponding diffusion matrix admits lower and upper bounds, uniform if 
< i? for any fixed R> 0. 

Step 2: Cut-of at a level |f |i = R. 

Let us introduce Markov time = inf{r > : |'w(t)|-^ > R}- Wc define 
the processes equal to for r e [0, tr] and satisfying the equation 

dv^ir) = dWr, T>fR, 

where Wr — ^ ^li ^ . Also, we define v^''^ to be equal to for r e [0, tr] 
and for r > tr satisfying the equation 

civ^'«(r) = U{yu{fR), n{fR))dWr , r > tr. . 

These processes have positive definite diffusion matrices uniformly in 7 and 
and 



38 



By Lemma 4.1 we have 

P{vl(r) ^ v^'^(r) for some < r < T} ^ 0, 
P{|vfc(r)| 7^ |vf (r)| for some < r < T} ^ 



as -R — cxD, uniformly in 7 and v. Therefore, it suffices to prove the lemma 
for Vfc replaced by with arbitrary R. ^ 

Step 3: limit 7^0. 

Denote a limiting (as 7 — >■ 0) law of v^'^ in C(0,T;M2) by and let 
Vfc(r) be any process such that its law equals CP. By construction, the 
relation holds £°{|vfc(T)| G Q} = C{\\^{t)\ e Q} for any Borel set Q C R. 
So it suffices to prove the lemma's assertion with replaced by Vjfc. 

The process v^' satisfies the relation 

r T 

^^l^r) = v«(0) + I i.,H(v^'^, V, s)ds + 1 5,,,^(v^'«, V, s)d(5l (7.7) 



with 



and 




-^k,R 



Bkj, s<tr, 



Bkj,R = \ ( I \ ^ , f 







where t/ = ?7(vfe(rR), v^(Ti?)). 

Denote in (7.7) the drift and martingale parts by ^^{t) and A4'^{r), 
respectively. Then 



r 







Distributions of the pairs {A"'{-), form a tight family of Borel mea- 

sures in C(0, T;R^). Consider a hmiting measure and represent it as the 



^Indeed, for any e > choosing first R so big that the probabihty in (7.6) is < e/2 and 
choosing next 6 = S{e) so small that the l.h.s. of (4.5), evaluated for replaced by Vj^, 
also is < e/2, we see that the l.h.s. of (4.5) is < e, if 5 is sufficiently small. 
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distribution of a process {A°{t),M^{t)). Then C{A^{-) + M^{-)} = so 
we can take for a process above the process Vfc(r) = A^{r) + A4^{r). Let 
Ti and T2 be arbitrary distinct point of [0,T] and Co - any positive number. 
The set G C(0,T;R2) : |v9(ri) - ip{T2)\ < Coin - ra]} is closed, thus 

]imsupP{\A^{n) - A^{t2)\ < Coin - T2I} 

7-0 (78) 
< P{|^°(ti) - ^V2)l < Colri - r2|} . 

Let us choose Cq = 2sup{\Ak{v)\ : \v\i < R}. Then 

\A^{n)-A^T2)\ < ^Co\n-T2\ + iy-'C{R)\{uAj) n [o,r]| . 

Since 

E|(UA,)n[0,T]| <p|sup|^(t)Ui >7~'|+E / XMr)\<2^dT, 

[o,r] Jo 

then it follows from (4.1) and Theorem 2.2.4 in [KrySO] that E|(UAj) fl 
[0,T]| — s> as 7 ^ 0. Therefore the limit in the l.h.s. of (7.8) equals 1, and 
we conclude that P{\A%ti) - A%T2)\ < Coin - T2\} = 1. That is, A^{t) is 

T 

Co-Lipschitz continuous and A^{t) = J B^{s) ds, where \Bq\ < Cq. 



We now turn to the martingale part. Since 

[0,T] 9 M^(t) e R^ 0<7<1, 

is a family of continuous square integrable martingales with respect to the 
natural filtration and uniformly bounded second moments, then the lim- 
iting process A4^{t) is a continuous square integrable martingale as well. 
Denote {{Ai'^))r the bracket (quadratic characteristics) of Ai'^. Accord- 
ing to Corollary Vl.6.7 in [JS87], {{M°))r = lim^_o((MT))^ . Since for 
V & {v : \v\hi < R} it holds 

with some Ci > 0, then the bracket {{Ai^)) satisfies the same estimate. In 
particular, d{{/A^))r = a{r)dr for some progressively measurable symmet- 
ric 2 X 2-matrix a(r) such that Ci Id < a(r) < cf^ Id, a.s. Then ~ 

T 

Jq a'^/'^{s) dM°{s) is a Wiener process in and M°{t) = J a^/^(s) dWg. 
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We have seen that for any u > and it! > 1 each weak hmit of the family 
v^'^(t) is an Ito process of the form 

r T 

Vfe(T)=Vfc(0) + J ds + J a'/^{s) dWs, 



where \B^{r)\ < Co and c]^^^Id< a^/^(r) < c^^^'^Id a.s., uniformly in t and 
I/. Since all the coefficients of this equation are uniformly bounded and 
the diffusion matrix is positive definite, the desired statement follows from 
Theorem 2.2.4 in [KrySO]. 



8 Appendix 

Here we prove the a-priori estimates, claimed in Section 1. 

Let F : H"^ — > R be a smooth functional (for some m > 0). Applying 
formally Ito's formula to F{u{t)), where u{t) is a solution, and taking the 
expectation we get 

^EF(«(t)) = E{VF{u),uu,, - V{u)) + ^uY.^lEd^^{u)[e,,e,]. 

s 

In particular, if F{u) is an integral of motion for the KdV equation, then 
(VF(m), V{u)) = and we have 

^ EF{u{t)) = uE{WF{u),u,,) + li^Yl blEd'F{u)[es, e,]. (8.1) 

s 

Since ||m||q is an integral of motion, then F{u) — exp((7||u||Q), < cr < |, 
also is an integral. We have: 

VF(«) = 2(Te^Wo^, d2F(M)[e,e] =2(7e^lHlo||e||2 + 4(7Vll"llo(M,e)^ 

So (8.1) implies that 

di 

where for r > we set 
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Denoting B — max 6^ and choosing a < {2B) ^ we get that 

AEe-^Wo < -i/(7E(e'^NIi(||^i||2_5^)^ < -i/aBo(Ee"Nlg _ 2e2'^^o). 

(JjV ^ / 

So the estimate (1.4) holds for all t > 0. In particular, for each > we 
have 

E\\u{t) 11^ <Mn^ Ca-^l'^ ■ { the r.h.s. of (1.4) ) . (8.2) 

The KdV equation has infinitely many integrals of motion Jm{u), m> 0, 
which can be written as 

m „ 

= \\U\\1 + Cr,mU^"''^ . . . U^"^^^ dx . (8.3) 

Q ™ ^ 



r=3 111 



Here the inner sum is taken over all integer r- vectors m = (mi, . . . m^) such 
that < rrij < m — 1 Vj and uii + . . .m^ = 4 + 2m — 2r (in particular, 
Jo = ||m||^). E.g., see [KP03], p. 209. 
Let us consider an integral as in (8.3), 

where / > 2, M > 1 and < m^- < - 1, O := //"^(M + //2 - 1) < 2 for 
some // > 2. Then, by Holder's inequality. 



|/|<M'"^^|l,,...M'"^')|l,,, p, = — <oo. 



M 

Applying next the Gigliardo-Nirenberg inequality we find that 

\I\<C\\ur,\\u\\i''. (8.4) 
Finally, evoking the Young inequality we get that 

\I\<S\\u\\l + Cs\\u\\l^ , V5>0. (8.5) 



We have 

m+2 



r=3 m 
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where mi H l-rrir — 6 + 2m — 2r. Due to (8.5) with S — 1/2, f — r and 

II — m + 1, 

h < -\ \\u\\l^, + C\\u\^P < -\ \\u\\l^, + C(l + \\u\\r^'^) . 
Next, 

r m 

Hence, 



r,m 



where f = r — 2 and mi + • • • + m,-. = 4 + 2m — 2r — mi — m2 =: M, M > 0. 
Note that |ej|c" = j" for each j and n. Assume first that r > 4 and M > 0. 
Then (8.4) imphes that 

h < 2f^ + ciix.||®+ir^+"^^iii.|ir'-® , 

with e = 2 - (3/2)'-+"»i+rn2 _ 

m+1 

By the Young inequahty, 

h < 2f^ + 5\\u\\i^i + csir'+'^'Wuro-'-^)^ 

It is easy to see that this estimate also holds for r = 4 and for M = 0. 

Using in (8.1) with F — the obtained bounds for Ii and I2 we get that 

I EUu) < -I ^nHL^i + ^i^E(i + ii^iir^^)) + ^ Yl 1^1""^' 

Choosing S — B^^ and using (8.2) we arrive at the estimate 
^ EJ„(m) < -i/E||m||^+i + uC^ , 
43 



where Cm depends on Bj^+i and M^m+i)- 

Applying (8.5) with ji — m to (8.3) we see that 

\ \\u\\l - C{1 + ||«||r) < Uu) < 2\\u\\l + C{1 + ll^lir). (8.6) 
Therefore 

where depends on the same quantities as Cm- We get that 

^JmHt)) < msix{EJmiu{0)),C'J 

for each t >0. Using (8.6) we obtain (1.5). 

Let us take any integers m > 0, > 1. By the interpolation inequality 
^ lklUfe||ii||o~'^- Therefore 

Using this inequality jointly with (8.2) and (1.5) we get the estimate (1.6). 
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